arXiv:1503.02458v2 [math.NA] 29 Aug 2016 


Shape Preserving Rational Cubic Spline Fractal Interpolation 


A. K. B. Chand and P. Viswanathan 
Department of Mathematies 
Indian Institute of Teehnology Madras 
Chennai - 600036, India 

Email: ehand@iitm.ae.in ; amritaviswa@gmail.eom 
Phone: 91-44-22574629 
Fax : 91-44-22574602 

Abstract 

Fractal Interpolation Functions (FIFs) developed through Iterated Function Systems (IFSs) 
offer more versatility than the classical interpolants. However, the application of FIFs in the 
domain of shape preserving interpolation is not fully addressed so far. Among various inter¬ 
polation techniques that are available in the classical numerical analysis, the rational interpo¬ 
lation schemes are well suited for the shape preservation and shape modification problems. 
Consequently, we introduce a new class of rational cubic spline FIFs that involve tension 
parameters as a common platform for the shape preserving interpolation and the fractal in¬ 
terpolation to work together. Suitable conditions on the parameters are developed so that 
the rational fractal interpolant retain the monotonicity and convexity properties inherent in 
the given data. With some suitable hypotheses on the original data generating function, the 
convergence analysis of the rational cubic spline FIF is carried out. Due to the presence of 
scaling factors in the rational cubic spline fractal interpolant, our approach generalizes the 
classical results on the shape preserving rational interpolation by Delbourgo and Gregory 
[SIAM J. Sci. Stat. Comput., 6 (1985), pp. 967-976]. Furthermore, for preserving shape 
of a data set wherein the variables representing the derivatives have varying irregularity, the 
present schemes outperform their classical counterparts. Several examples are supplied to 
support the practical value of the method. 

This paper was submitted to SIAM Journal of Numerical Analysis on April 14, 2013 
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1 Introduction 

In the classical Numerical Analysis, there are several interpolation methods that can be applied to 
a specific data set, according to the assumptions that underlie the model we investigate. However, 
if the given data set is more complex and irregular (for instance, data sampled from real-world 
signals such as financial series, seismic data, and bioelectric recordings), then the traditional in¬ 
terpolants may not provide satisfactory results. To address this issue, Barnsley [[D introduced 
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a class of functions called FIFs using the notion of IFS. FIFs aim mainly at data which present 
details at different scales or some degree of self-similarity. These characteristics imply an irreg¬ 
ular structure to the interpolant. The main differences of FIFs from the traditional interpolants 
reside in: a) the definition in terms of functional equation which implies a self similarity in small 
scales; b) the iterative construction instead of using an analytical formula and c) the usage of 
some parameters, which are usually called scaling factors that are strongly related to the fractal 
dimension of the interpolant. Later, Barnsley and Harrington ||4l observed that if the problem is 
of differential type, then the parameters of the IFS may be chosen suitably so that the correspond¬ 
ing FIF is smooth. This observation initiated a striking relationship between the classical splines 
and the fractal functions. Smooth FIFs (fractal splines) constitute an advance in the techniques 
of approximation, since the classical methods of real-data interpolation can be generalized by 
means of smooth fractal techniques (see, for instance, [r71llO[l34[l35H l. Further, if the experimental 
data are approximated by a C^-FIF S, then the fractal dimension of provides a quantitative 
parameter for the analysis of the data, allowing to compare and discriminate experimental pro¬ 
cesses. Though FIFs are primarily applied for a self-affine/self-similar data set, their extensions 
namely hidden variable FIFs [O and coalescence hidden variable FIFs (HI can be used to simu¬ 
late curves that are non-self-affine or partly self-affine and partly non-self-affine in nature. Due 
to these versatility and flexibility, the theory of FIFs has evolved beyond its mathematical frame¬ 
work and has become a powerful and useful tool in the applied sciences as well as engineering. 

The central focus of interpolation (traditional or fractal) is to construct a continuous function 
that fits the given points obtained by sampling or experimentation. However, to obtain a valid 
physical interpretation of the underlying process, it is important to develop interpolation schemes 
that inherit certain properties from the prescribed data set. Examples of few such prevalent fea¬ 
tures are positivity, monotonicity, and convexity. Constraining the range of an interpolation func¬ 
tion so as to yield a credible visualization of the data by adhering to these intrinsic characteristics 
is generally referred to as shape preserving interpolation. There are multitudes of classical in¬ 
terpolation methods that honour shape properties inherent in the data. In what follows, we shall 
provide some pioneering works in this field. 

Research on shape preserving interpolation has been originated with some existence-type 
results by Wolibner (421 and Kammerer (291 . These results do not provide any additional in¬ 
formation on the shape preserving polynomial. A constructive approach to the shape preserving 
interpolation using hyperbolic tension splines was popularized by Schweikert (38l . Main issues 
connected with the hyperbolic tension splines are: (i) development of an automatic algorithm for 
the choice of free parameters is complicated; (ii) it is computationally complex to work with, 
especially for very large/small values of tension parameters involved in it. Polynomial splines 
gain shape properties either (i) by addition of extra knots (see (211 [33l 13^1401 ). which may not 
be effective in terms of computational economy, or (ii) by perturbing derivative parameters (see, 
for instance, (2211231 ). which make the method unsuitable for Hermite data, where the given 
derivative values are also to be interpolated. In shape preservation and shape control, rational 
splines provide an acceptable alternative to the polynomial/hyperbolic splines. Wide applica¬ 
bility of the rational interpolants may be attributed to their ability to receive free parameters 
(which may be used for shape control) in their structure, ability to accommodate a wider range 
of shapes than the polynomial family, excellent asymptotic properties, capability to model com- 
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plicated structures, better interpolation properties, and excellent extrapolating powers. Gregory 
and Delbourgo popularized the shape preserving rational interpolation methods through a series 
of papers [I16II18II191I24II25H . These works stimulated a large amount of research in the direction of 
shape preserving rational spline interpolation. For brevity, the reader is referred to II131I261I271|J71 . 

These non-recursive shape preserving interpolation techniques, in general, produce smooth 
interpolants whose derivatives are also smooth except possibly at some finite number of points. 
However, in practice, there are many situations where the variable in the data possesses certain 
shape properties, and at the same time, variables representing the derivatives may be irregular. 
For instance, a sphere falling in a wormlike micellar solution does not approach a steady terminal 
velocity, instead, it undergoes continual oscillations as it falls f[2^ . Hence, to simulate the dis¬ 
placement and velocity profiles of such motions, monotonicity/positive interpolants with varying 
irregularity (fractality) in the derivatives may be advantageous. Similarly, in nonlinear control 
systems iHTIl (say, for instance, the motion of a pendulum on a cart) monotonicity/convexity pre¬ 
serving interpolants with varying irregularity in the second derivative may be desirable for the 
study of acceleration. Therefore, it is useful to develop smooth FIFs (which are known to have 
fractality in their derivatives) that retain the intrinsic shape properties of the data set. From the 
knowledge gained from the classical shape preserving polynomial interpolation techniques, it 
is felt that preserving fundamental shape properties via polynomial FIFs would be difficult or 
impossible. Thus, for an initial exposition of FIFs to the field of shape preserving interpolation, 
rational FIFs seem to be an appropriate medium. 

With these motivations, the capability of FIFs to generalize smooth classical interpolants, 
and the effectiveness of rational function models in shape preservation are intertwined to pro¬ 
vide a new solution to the shape preserving interpolation problem from a fractal perspective. 
We construct a -rational cubic spline FIF with one family of shape parameters in section [3Tj 
To demonstrate the effectiveness of the rational cubic spline fractal interpolation scheme, the 
convergence results are discussed in section 13.21 The rationale behind selecting a rational FIF 
with shape parameters instead of widely studied polynomial FIFs is the following. If the scaling 
factors tend to zero and the shape parameters tend to infinity, then the rational FIF converges to 
the piecewise linear interpolant for the data. This tension effect ensures that the FIF can be used 
for shape preserving interpolation. Section |4] provides an automatic selection of the parameters 
that culminate in interactive algorithms to preserve monotonicity/convexity of the data. These 
algorithms take full advantage of the flexibility that the fractal splines permit. By suitable choice 
of the shape parameters that verify the monotonicity condition, our cubic spline FIF reduces to a 
lower degree rational spline FIF that generalizes the classical rational spline interpolant studied 
in [[25l . With special choices of the shape parameters satisfying convexity condition, our cu¬ 
bic/quadratic rational FIF reduces to lower order form, which provides the fractal generalization 
of the classical rational interpolant discussed in ifT^ . Again, by proper choices of the scaling 
factors and the shape parameters, our rational cubic FIF degenerates to the classical rational cu¬ 
bic interpolating function introduced in lfT9l . Therefore, the present paper offers a novel idea of 
setting a common platform for the fractal interpolation and the shape preserving interpolation 
to operate together, and in the process collectively generalizes three different classical rational 
interpolation schemes available in the literature. In section [51 some remarks and possible ex¬ 
tensions are made. The effectiveness of our shape preserving fractal interpolation schemes is 
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illustrated with suitably chosen numerical examples and graphs in section 0 


2 Basics of Polynomial Spline FIF 

To equip ourselves with the requisite general material for the construction of the desired rational 
spline FIF, we shall reintroduce the polynomial spline FIF. A complete and rigorous treatment 
can be found in 

2.1 Fractal Interpolation Functions 

Let {(xj, Di) e / X M : i = 1, 2,..., A^} be a real data set, where xi < X 2 < • • • < xat is a 
partition of / = [xi,X 7 v]- Set, K = I x D, where D is a large enough compactum in M. Let 
J = {1, 2,..., — 1}, and Li : I —)■ Jj = [x,, Xj+i] be affine maps satisfying: 

Li(xi) = Xi, Li(xAr) = Xi+i, i e J, (2.1) 

and Fi : K —)■ D be continuous functions such that: 

Fi{xi,yi) = Vi, Fi{xN,yN) = 2/i+il 
\Fi{x,y)-Fi{x,y*)\<\ai\\y-y*\j' 

where (x,j/), {x,y*) G iL, 0 < \ai\ < k < 1 for alH e J, and k is a fixed real constant. 
Define Wi{x,y) = [Li{x), Fi{x,y)) for z e J. It is known [HI that there exists a metric on 
equivalent to the Euclidean metric, with respect to which Wi, i G J, are contractions. The 
collection {iL; w*, z e J} is termed as an Iterated Function System (IFS). On 'H(iL), the set of 
all nonempty compact subsets of K, endowed with the Hausdorff metric, define a set valued map 
W (A) = [J Wi{A). Then, VF is a contraction map on the complete metric space 'H{K). Thanks 

i£j 

to Banach Fixed Point Theorem, there exists a unique set G E 'H{K) such that W{G) = G. The 
set G is termed as the attractor or deterministic fractal corresponding to the IFS {iL; zxj, z G J}. 
The definition of a FIF originates from the following proposition: 

Proposition 2.1. (Barnsley /[T]/) The IFS {K;Wi, i E J} has a unique attractor G such that G 
is the graph of a continuous function /:/—)■ M which interpolates the data {(xj,z/j) : z = 
1,2,..., A^}, i.e., G = {(x,/(x)) : x G /} and f{xi) = Pifori = 1,2,..., A^. 

The function / in Proposition 12.1! is called a fractal interpolation function corresponding to 
the IFS {K]WiG ^ J}- The adjective/ractoZ is used to emphasize that G = graph(/) may 
have noninteger Hausdorff-Besicovitch dimension. But / may be many times differentiable (see 
section |2^ . Since G is a union of transformed copies of itself, i.e., G = [J wfG), an alternative 

i€J 

name for a fractal function could be a self-referential function. The characterization of a graph of 
a FIF by an IFS leads to a recursive construction of / using the following functional equation [HI : 

f(x) = F,{L-\x)JoLf\x)), xEf, lEJ. (2.3) 
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The following special class of IPS has received wide attention in the literature: 

Lj(x) = ttiX + bi 

Fi{x,y) = aty + Ri{x) 

where i G J, are parameters satisfying \ai\ < k < 1 and : / —)■ M, i G J, are suitable 
polynomials satisfying (12.21) . The multiplier a* is called a scaling factor of the transformation 
Wi and a = (ai, 02 ,..., oat-i) is the scale vector of the IPS. As mentioned earlier in the intro¬ 
ductory section, a PIP is defined in a constructive way through iterations instead of descriptive 
one, usually a formula, provided by the classical methods. Consequently, evaluation of a PIP 
at a given point needs, in general, an iteration process. An explicit representation (in terms of 
an infinite series) of a PIP / corresponding to a general IPS (12.41) defined on J = [0,1] is given 
in m. Por a representation of / as the uniform limit of a sequence of operators, the reader is 
referred to flUl. 


G J, 


(2.4) 


2.2 Polynomial Spline FIFs 


Por a prescribed data set, a polynomial PIP with C''-continuity is obtained as the fixed point 
of IPS (12.41) . where the scaling factors a, and the polynomials Ri involved in (12.41) are chosen 
according to the following proposition. 


Proposition 2.2. (Barnsley et al. /Ej/) Bet x\ < X 2 < ■ ■ ■ < xn and Li(x) = atx + bi, i ^ J, 
be the affine functions satisfying (I2.il) . Let Ffx, y) = aty + Ri(x),i G J, satisfy (12.21) . Suppose 
that for some integer p >0, |Q!j| < of and Ri G C^{I), i E J. Let 


Fi,k{x,y) = 


ttiy + Rf^ (x) 


_ K-iM 

■) yi,k u ) yN,k u 1 

al — ai a.'ff-i ~ <xn-i 




k = 1,2,...,p. 


IfFi_i^k{xN, yN,k) = Fi^k{xi,yi^k) for z = 2, 3,..., - 1 and k = 1, 2,... ,p, then the IFS 

(^Lfx), Fi{x, y )), z G J} determines a FIF f G CP[xi,X]\f], and is the FIF determined 
by {Li(x),Fi^k{,x,y)),i e J] for k = 1,2,... ,p. 


3 Rational Cubic Spline FIFs Involving One Family of Shape 
Parameters 

In section [Xn we construct (Irrational spline PIPs where the inhomogeneous terms are rational 
functions with cubic numerators and preassigned quadratic denominators. In section 1X21 error 
analysis of the rational cubic spline PIP is given with the assumption that the data defining 
function / G C^. Purther, by admitting a relatively weaker condition on the data generating 
function /, namely / G C^, the uniform convergence of the classical rational cubic spline is 
established. This serves as an addendum to the convergence results by Delbourgo and Gregory 
IIT 9 I . and it is utilized to establish the uniform convergence of the developed rational cubic spline 
PIP. 
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3.1 Existence of C^-rational Cubic Spline FIF 


Theorem 3.1. Suppose a data set {{xi, pi, di) : i = 1, 2,, A^} is given, where xi < X 2 < 
■ ■ ■ < xjsi. Consider the rational IPS {{Li{x),Fi{x,y) : i G J} where Li{x) = OiX + bi and 
Fi{x, y) = Uip + Ri{x), Ittil < KOi, 0 < K < 1, i G J. Further, let Ri{x) = where Pi{x) is 
a cubic polynomial and Qi{x) 0 is a preassigned quadratic polynomial such that Fi{xi, yi) = 

Hi, Fi{xN,yN) = Pi+i are satisfied. With Fi^i{x,y) = i^t Fi^i{xi,di) = di and 

FfixN, (^iv) = (ii+i. Then a C^-rational cubic spline FIF S satisfying S{xi) = pi, SR'>{xi) = di, 
i = 1, 2,..., iV exists, and it is unique for a fixed choice of the shape parameters and the scaling 
factors. 

Proof. Set I = [xi,xn], and R = Xi+i — Xi, i E J. Consider the IPS T = {{Lfx), Ffx, y)) ■ 
i G J} where Lfx) = OjX + bi satisfy (12.11) . and Ffx^y) = aip + Rfix) fulfill the join-up 
conditions Ffixi^yf = Pi, FfxN.PN) = Vi+i- Let 


Rfx) = R* [xi + 9{xn - a:i)) 


Afi - efi + Bie{i - efi + -e) + DiO^ 

1 + (r, - 3)0(1 - 9) 


Consider R := {g E C{I)\ g{xi) = yi and g{x]\f) = yj^}. The uniform metric d{g,h) : = 
sup{|fi'(a;) — h{x)\ : x E 1} completes R. The IPS Xinduces a contraction map T : R ^ R 


{Tg){Li{x)) := Fi{x,g{x)), x El. 


The contraction map T has a unique fixed point S E R, which satisfies the functional equation: 

S(Li{x)) = Fi{x,S{x)), 

^, A(1 - 9f + Bfi{l - 9y + C^9\l - 0) + (3.1) 

= +- l + (r.-mi-0) -■ 

The conditions Ffxi, yf = Pi, Ffix^, y^) = Pi+i can be reformulated as the interpolation con¬ 
ditions S{x.i) = Pi, S'(a;j+i) = Pi+i, i E J. Note that: (i) the affine map R satisfies Lfixi) = Xi 
and L^xn) = Xi+i, (ii) x = xi and x = xn correspond to 0 = 0 and 9 = 1 respectively. There¬ 
fore, the interpolatory conditions determine the coefficients Ai and Di as follows. Substituting 
X = xi'm (13.11) we get 


S{Li{xi)) = aiS{xi) + Ai, 


Vi — OtiPi -f Ai, 

Ai = Pi- aipi. 


Similarly, taking x = x^ in (13.11) we obtain Di = pi^i — Uip^. 

Now we make S' G by imposing the conditions prescribed in Bamsley-Harrington theorem 
(see Proposition l2.2l) . 

Assume \ai\ < KOi, i E J, where 0 < k < 1. We have Ri E Adhering to the notations of 

Propositionwe let: 


RAx,y) 

2 / 1,1 = di, pn^i 


aiP + Rf\x) 

Qji 

dN, FiRxi,di) 


di) Fi i(^x]\f, d]\f^ 


dj+i; i E J. 
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Then, by Proposition l2.2l the FIF S belongs to the class C^(/). Further, is the fractal function 
determined by the IFS X* = {[Li{x), Fi^i{x,y)) : i G J}. Consider X”* := {g G C{I) : g{xi) = 
di and g{xN) = dN} endowed with the uniform metric. The IFS X* induces a contraction map 
T* : F* -)■ F* 

{T*g*){L,{x)) ■.= F,,,{x,g*{x)), xgX 

The fixed point of T* is Consequently, satisfies the functional equation: 




(3.2) 


The conditions on the map X i, namely X i(a;i, di) = di and Fi^i{xN, d]sf) = X+i can be refor¬ 
mulated as the derivative conditions = X and = X+i, i G J. 

Applying x = xi in (13.21) we obtain 


Bi 


TiAi 


S^^\Li{xi))ai = + 

Xat — Xi 

diai{xN - xi) = aidi{xN - xi) + X - ri{yi - aiyi), 
^ Bi = [viyi + hidi] - ai[riyi + di{xN - a;i)]. 

Similarly, the substitution x = in (13.21) yields 

Ci = [viyi+i - hidi+i] - ai[riyN - dN^XN - a;i)]. 

Therefore, the desired C^-rational cubic spline FIF is described as: 

Fix) 


S{Li{x)) = aiS{x) + 


Qiix)' 


(3.3) 


Pi{x) = P*{e) = {yi - aiyi){l - 6)^ -f {[rtyi + hidi] - ai[riyi + di{xN - a:i)]}6'(l - 9)^ 
+{[riyi+i - hidi+i] - airiy^ - dN{xN - Xi)]}X(l - 9) + (j/^+i - aiyN)9^, 

Q,{x) ^ Q*{9) = l + {n- 3)0(1 -9), 9 = 

The parameters r* > — 1 can be effectively utilized for the shape modification and shape preser¬ 
vation of the C^-rational cubic spline FIF, and hence referred to as the shape parameters. 

Since the FIF S in (13.31) is derived as a solution of the fixed point equation Tg = g, the solution 
is unique for a given choice of the scaling factors and the shape parameters. □ 

Remark 3.1. Using the notations 

0(1 - 0)X{(20 - l)Ai + (1 - 9)di - 0X+i} 


X:= 

Fi := 


1 + (a ~ 3)0(1 — 0) 

0(1 - 0){(20 - l)(?/jv - yi) + (1 - 0)(xjv - xi)di - 0(xjv - xQdjy} 

1 + (r^ - 3)0(1 - 0) 
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the rational cubic spline FIF in 0.3D is rewritten as 


S{Li{x)) = aiS{x) + {(1 - e)yi + %+i + Ei} - ai{{l - 9)yi + OyN + Fi] (3.4) 

As the shape parameters ri are increased, Ei and E^ converges to zero. Thus from 0.4D . it follows 
that as the scaling factors a* —)• 0 and the shape parameters rj —)■ +cxd, the rational cubic spline 
FIF S converges to the piecewise linear interpolant corresponding to the given data set. This 
tension effect ensures that the proposed rational cubic spline FIF can be used to construct shape 
preserving interpolants. 

Remark 3.2. If ai = 0 /or all i E J, then S reduces to the piecewise defined -rational cubic 
spline s discussed in H^l . Therefore, S can be considered as an extension of the powerful 
rational cubic spline interpolant. To illustrate this we proceed as follows. With a* = 0 for all 
i E J, 0.31) reduces to 


S{Lfix)) = 


yfl - Of + {rtyt + hidfiOf - Of + {riPi+x - hidi+x)0‘^{l - 6) + Pi+iO^ 

1 + (r, - 3)0(1 - 0) 


i E J. 
(3.5) 


Since 


Xjsi — Xl 


from ^3.51) . for x E f = [xi, Xj+i], we have 


S{x) = s/x) 


yfl - fif + {riPi + hidfjfiil - fif + {viPi+i - - /) + yi+i(f 


1 + (r, - 3)/(l - /) 
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(3.6) 


where fi = is a localized variable. The rational cubic spline s E C^{I) is defined by 
s| = Si, i E J. With Ti = 3 and the scaling factors satisfying laf < noi, W i E J, our 
discussion on the rational cubic spline FIF gives a simple constructive approach to the -cubic 
Hermite FIF. Again, when = 0 and = 3,'i i E J, the proposed rational cubic spline FIF 
recovers the classical piecewise cubic Hermite interpolant. 


3.2 Convergence Analysis of Rational Cubic Spline FIFs 

With mild conditions on the scaling factors, we establish that the rational cubic spline FIF S 
possesses the same convergence properties as that of its classical counterpart s. Since S does 
not possess a closed form expression, standard methods such as Taylor series analysis, Cauchy 
remainder form, and Peano Kernel theorem cannot be employed to establish its convergence. In¬ 
stead, we derive the convergence of S to the original function / using the convergence results for 
its classical counterpart s and the uniform distance between S and s via the triangle inequality: 

||/-^||oo<||/-s|U + ||s-5|U. (3.7) 

Theorem 3.2. Let S and s, respectively, be the rational cubic spline FIF and the classical ra¬ 
tional cubic spline for the original function f E C‘^(/) with respect to the interpolation data 
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{{xi, Di) : i = 1, 2,..., N}, and let S'^^\xi) = = di. Suppose that the rational function 

Ri involved in the IPS generating the FIF S satisfies | | < Z^for \Ti\ G ( 0 , nai), all i & J, 

and for some real constant Zq. Then, 


h 


11/ - 5||oo <£; max{|j/f^ - di\, - di+f} + + jmaxiri - 3|) 


.( 1 ) 


'Ac i€J 


384c 


4 max|ri - 3|(/i^||/(^)||oo + 3/i^||/^^^|loo) } + 1^^122^ 
ieJ i — a 


4 idJ 

Zo) 


(3.8) 


where |a|oo = max{Q;j : i G J}, h = max{/ij : i G J}, and c = min{cj : i G J} wiY/i 


3^, -1 < < 3, 

1, Vi > 3. 


(3.9) 


Proof For a prescribed set of data and a* satisfying |aj| < nOi, i E J, the rational cubic spline 
FIF S G C^(J) is the fixed point of the Read-Bajraktarevic operator Tq,: 


{TaS){x) = aiS{L. ^(x)) + Rfai, 0), 


(3.10) 


p:iai,4>) 


— X G [xi, Xi+i], i E J, with P* and Q* as in 


where Rfa^, f) = ^ 

(13.31) . Note that the subscript a is used to emphasize the dependence of the map T on the scale 
vector a. The coefficients of the rational function Ri depend on the scaling factor a*, and hence 
Ri can be thought of as a function of and f. The interpolants S and s are fixed points of 
with a 7^ 0 and a = 0 respectively. We know ifTOl that for x E [x*, Xj+i], 


I/M - s(i)| <^max||!/‘'' - iil, l!/,™ - (ij+i 

+ 4|n-3|(A||l/W||+3ft?||/W||)}, 


384c,- 


ht\ma + 


fi - 3| 


(3.11) 


where = f^^\xi) for i = 1 , 2 ,..., A^, and ||.|| denotes the uniform norm on [x,, Xi+i]. 
For a fixed choice of scale vector a 7^ 0 and for x G [x,, x,+i], from (13.101) we obtain: 


|TaS'(x) - Tas{x)\ = {aiS{Li ^(x)) + Ri{ai,(j))} - {a/s(Lj ^(x)) + Ri{ai,(j))} 
^ 10:10011*5* s||oo. 


From the above inequality we deduce: 


llTaS* - T„s|loo < 1010011 *5* - s||oo. 

Let X G [xj, Xj+i] and a 7 ^ 0. Using (13.101) and the Mean Value Theorem: 

|r„s(x) - Tos(x)| = {ais{L~^{x)) + Ri{ai,(j))} - Ri{0,(j)) , 

dRi{Ti,(j)) 




+ Oj 


< |o,; 


00 

s||oo + Zq), 


dai 


(3.12) 
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Thus, 


ooV oo 


||T„s - Toslloo < |a 

Using (13.121) and (13.131) . 

11*5' - s||oo = ||7'„S' - TqsIIoo < \\TaS - T„s||oo + \\TaS - ToS||oo, 

— |*^|oo||*5' slloo T |®|oo(||'S|loo T ^0 )) 


(3.13) 


which simplifies to 


11*5 - s||oo < 


|q^|oo(||'5||oo t 

1 I I no 


(3.14) 


The required assertion follows from (13.71) . (13.111) . and (13.141) . However, in what follows, we find 
an upper bound for ||s||oo and estimate Zq, if not optimally, at least praetieally. 

Let us introduee the notations: \y\oo = max{|?/j| : 1 < z < N}, \d\oo = max{\di\ I < i < N}, 
and |r|oo = max{|rj| : i G J}. From (13.61) . for x G [xi,xn]. 


|s(a:)| < 


max{|P"((/))| : i e J,0 < (f) < 1} 


min{|( 5 *( 0 )| : z G J, 0 < 0 < 1 } 
where Pi*{4>) is the numerator in (13.61) . Using extremum ealeulations of polynomials, 

\Pi*i<P)\ < |yi|(l - 4’f + {\ri\\yi\ + hi\di\)(j){l - (j)f + (\ri\\yi+i\ + hi\di+i\)(p‘^{l -</>) + \yi+i\(p^, 
=)> max^jP**< max{|yj|, |yi+i|} + ^(k*l max{|yi|, |yj+i|} + himax{\di\, |di+i|}^, 

^ max |F-**((())| < I 2 /I 00 + 7 (|r|oo| 2 /|oo +/i|d|oo), 
ieJ, 06 ( 0 , 1 ] 4 

and \Q*{(j)) \ = Qi(0) > U. Therefore, 

|?/|oo T 4 ( I'^loo |?/|oo T ^M|oo) 


(3.15) 


II "5II 00 ^ . r . J., 

minjcj ; z G J| 

Now, from (13.31) and (13.101) . for x G [xj, Xj+i], 

dR, _ P,(0) 
dai Qi{(j))' 

where Pi{(j)) = - (j)f + {rtyi + {xn - xi)di)(j){l - (j)f + {rtyN - dN{xN - - 

0) + Using similar extremum ealeulations. 


^R^ 

/ ^ max{|z/i 

) 2/Af }(1 + i 

r 

00 ) + i 

/ max{ di 

, Miv} 

dai 


min{ 

Cl 

: z G J} 



Vz G J, 


where |/| = xn — xi. 


(3.16) 

□ 
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Due to the principle of construction of a smooth FIF, for S to be in the class we impose 

I a, I < KQi = ■ Hence, |a|oo < x^-xi ’ consequently S converges uniformly to the 

original function when the norm of the partition tends to zero. The following convergence results 
are direct consequences of Theorem 13■2[ 

Corollary 3.1. Let S be the rational FIF with respect to the data points {(xj, |/j) : i = 1, 2,..., N} 
corresponding to the original function f G C^(/). Suppose dj, i = 1, 2,..., A^, and Vi, ai,i ^ J 
are chosen accordingly. 

(i) Ifyf'’ - di = 0{hi) = - di+i, r* > -1 and \ai\ < of, then ||/ - S'||oo = 0{li^). 

(ii) Ifyf'‘ -di = = yf^^-di+i, rj-3 = 0{hi) and\ai\ < af then ||/-S'||oo = 0{h^). 

(in) Ifyf^ -di = 0(h^) = -d^+i, = 0(hj) and\ai\ < af, then |l/-S'||oo = 0{h^). 

The above theorem and corollary show that r* should ideally be such that ri — 3 = 0(h‘f). 
Later we shall consider how r* can be chosen to preserve the data monotonicity/convexity, whilst 
maintaining this optimal requirement. 

Following the convergence results for the classical rational cubic spline s studied by Gregory and 
Delbourgo, we assumed that the data generating function / is in class C"^(/). Now we establish 
the uniform convergence of s with a weaker assumption / G C^{I), and use it to deduce the 
uniform convergence of rational cubic spline FIF S as in the previous case. 


Theorem 3.3. Let S and s, respectively, be the rational cubic spline FIF and the classical ratio¬ 
nal cubic spline interpolant for the original function f G C^{I) with respect to the interpolation 
data {(xi,yi) : i = 1,2,... ,N}. Suppose that the rational function Ri involved in the IFS 
generating the FIF S satisfies [ \ < Zq, |rj| G (D, nafi, for all i G J, and for some real 

constant Zq. Then, 


ll/-^l|oo< 


— h\d\oQ + —uj{f] h){\r\oQ + 4) + 


a 

oo(||s| 

1 oo 

+ Zq) 

1 - 

a 

OO 


where c = min{cj : i G J}, and uj(f; h) := sup {\f{x) — f{x*)\ : x, x* G /} is the modulus 

\x—x*\<h 

of continuity of f. In particular, S converges uniformly to f ^ ihe mesh norm tends to 

zero. 


Proof. Observe that 

Qfe) = i + {r,-3)e{i-e), 

= (1 - 9)^ + rfiil - of + ne\l -9) + 9\ 
For X G [xi,Xi+i], 

P*(9) 


fix) - six) =/(x) - 


Q^i9) 

1 


(1 - ^f{fix) - Vi) + rfiil - 9)^{fix) - yi) + 


1 + (r, - 3)0(1 - 9) L 

- 0)(/(x) - yt+i) + 0^(/(x) - yi+i) - hfi^il - 0)^ + hidi+i9‘^il - 9) 
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Therefore, local error of the interpolation is given by 


\f{x)-s{x)\ <-< \ f{x) -|/i|[(l-0)^ + \ri\e{l-ef] + |/(a;) - yi+i\[\ri\e\l - 6) + 6^] 

Ci ^ 

+ hi\\di\6{l — OY + — 9 


<- 


d L 4 


^Vil + l)u{f]h) + ^/ijmax{|di|, \di+i\} 


Consequently, we have the following uniform error bound for the classical rational spline s: 

(3.17) 

Now (13.71) coupled with (13.141) and (13.171) proves the theorem. □ 


11/ - s||oo < ^h\d\oo + h){\r\^ + 4). 


4 Shape Preserving Rational Fractal Interpolation 


4.1 Monotonic Data 

For the sake of simplicity, let us assume that the given data set is monotonically increasing, i.e., 
2/1 < 2/2 < ■ • • < l/AT, and consequently Aj = > 0 V i G J. For a monotonic increasing 

interpolant S G C^, it is necessary that the derivative parameters satisfy dj > 0, i = 1, 2,..., iV. 
From elementary calculus, we know that a differentiable function S is monotonic increasing on 
/ if and only if > 0 for all x G /. Calculation of S^^\Li{x)) from (13.31) and further 

simplification give: for x G /, 


di 


- 0) + Vi9\\ - ey + v)0(i 

[1 + (r, - 3)0(1 - 9)y 


I ' » 7 


(4.1) 


Q/ • 

where = di+i - ^(x^v - xi)dAr, 

hi 

2ai 

Si = 2{riAi - di) - -j^[ri{yN - 2/i) - di{xN - aii)], 

cx ■ 

Ui = (Y + 3)Ai - ri{di + di+i) - T^[(r,^ + 3)(2/Ar - 2 / 1 ) - ri(xAr - Xi)(di + d^v)], 


hi 


2cyi 


Vi = 2{riAi - di+i) - ^[ri(2/jv - 2/i) - d^ixM - a:i)], 


di 


Wi = di - -r{xN - xi)di. 
hi 

Since the rational cubic spline FIF is defined implicitly and recursively, to maintain the positivity 
of in the successive iterations and to keep the desired data dependent monotonicity condition 
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to be simple enough, we assume a* > 0 for all i e J. Our predilection to the nonnegativity as¬ 
sumption on the scaling factors is attributable to reasons of convenience rather than of necessity. 
Then, for i G J and an arbitrary knot point Xj, sufficient conditions for [Li{xj)) > 0 are: 

Ti > 0, > 0, Ui >0,Vi> 0, Wi > 0, (4.2) 


where the necessary condition on the derivative parameters are assumed. 

Now Ti > 0 — xi)dN < di+i. Observe that if dN = 0, then T* > 0, i G J, follow 

directly from the assumption on the derivative parameters. Otherwise, we impose the following 
condition on the scaling factors: 


a,: < 


di-\-ihi 


dnixN -xi)' 


(4.3) 


Similarly IT* > 0, whenever the scaling factor satisfies 


Oii < 


dihi 


di{xN -xi)' 


(4.4) 


Again 


S'i > 0 hi{riAi - di) > ai[ri{yM - Vi) - di{xN - a:i)], 

1/. > 0 hi{ri^i - di+i) > ai[ri{yN - yi) - dN^XN - xi)]. 


If it is assumed that 


ri[hiAi - ai{yN - yi)] > hi{di + di+i) - ai{xN - xi){di + ^ at ), (4.6) 

then it follows from (I4.3I) - (I4.4I) that 

ri[hiAi - ai{yN - yi)] > hdi - aidi{xN - xi) > 0, 
ri[hiAi - ai{yN - yi)] > Kdi+i - aid^ixM - Xi) > 0 . 

In view of (14.51) . the above inequalities imply Si > 0 and Vi > 0. 

Assume that hiAi — ai{yN — yi) > 0, i.e.. 


_ 

Vn-Vi 

From (14.61) . we have 


r‘^[hiAi - ai{yN - yi)] > rihi{di + di+i) - riai{xN - Xi){di + d^). 


(4.8) 


From (14.81) and S/ijAj > 3ai{yN — |/i), it is easy to verify that Ui > 0. From (14.31) . (14.41) and 
(14.71) . for a monotonicity preserving rational FIF, it suffices to choose a*, z G J, according to: 


0 < Oj < min 


dihi 


Ktti 


di-\-\hi 


Aihi 


di{xN - Xi) ’ d^ixN - Xi) ’ Z/AT 


Vi 


(4.9) 


13 
















Here the first term in the braees arises due to the C^-eontinuity of S. After ehoosing the sealing 
faetor ai aeeording to (14.91) . the shape parameters is seleeted to fulfill inequality (14.61) . and 
these eonditions are suffieient for [Li{xj)) > 0. As [a;i,XAr] is the attraetor of the IPS 
{M; Li{x),i G J}, by the reeursive nature of the rational fraetal funetion, (Li(a;j)) > 0 for 
alH e J and for every knot point Xj imply that (a:) > 0 for all x G /. 

With the neeessary eondition d* < 0, i G J, assumed to hold, an analogous proeedure applies for 
a monotonie deereasing data set. 

Remark 4.1. A* = 0, then we take «* = 0 for the monotonicity of the rational cubic spline 
FIF. Also in this case, di = dj+i = 0. Consequently, S^Lfx)) = Pi = Pi+i ,i.e., to say that S 
reduces to a constant on the interval [xi, Xi+f. 

Remark 4.2. When all ai = 0, the rational cubic FIF reduces to the classical rational cubic 
spline s. In this case, condition ( 14.91) is obviously true, and the condition (14.(51) reduces to 


A > {di + di+i)//Ai, i G J. 


(4.10) 


VN-yi' 


Thus ( l4.iOD is a sufficient condition for the monotonicity of the classical rational cubic spline s 

[ mi, P- 970]. 

Remark 4.3. For a given strictly monotonic data, we select ai satisfying ( 14.91) with ai < 
and then fix the shape parameters according to 

_ ^ hi{di + di+i) - afxN - xi) ( 4 i + dAr) 

h^Ai-afivN-yi) ’ ^ 

so that the monotonicity condition (14.(51) is satisfied. With these choices of the IFS parameters, 
the rational cubic spline FIF reduces to the rational quadratic FIF: 

P^{d) 


S{Li{x)) = aiS{x) + 


Q^{ey 


(4.12) 


where Pi{e) = {yi+i-aiyN)Aie‘^+l3ii^{yidi+i+yi+idi)ai-ai[yi+idi+yidN+ai{yNdi+yidi+i)\ 

+af{yNdi + ?/i4jv)|6'(l - 9) + {pi - aipfAfl - Of, 

Qi{d) = Aj9^ + fii 'yafidi + dj+i) — afidi + c/at)] 9(1 — 9) + Aj(l — 9)^, and 

fii = --r. For Ai = 0, we choose ai = 0, and define S{Li{x)) = Pi = 

{Pi+i - Pi) - afip N - V i) 

Pi+i. If all a* = 0 in (14.721) . then the corresponding rational quadratic FIF reduces to the 
classical monotonic rational quadratic interpolant studied in detail in / l25l/ . In this degenerated 
case, the necessary condition 4* > 0 is also sufficient for the monotonicity of the interpolant 
(see mi)- 

Remark 4.4. For the shape parameters specified in (14.771) and |ai| < af we get ri — 3 = 0{hl). 
Consequently, from corollary 3.1 it follows that A4.11\) is a good choice of Vi since the optimal 
0{hy bound on the interpolation error can be achieved. Hence, for a monotonic rational cubic 

spline FIF with an optimal error bound, we choose 0 < < min , /t* 

K, n* G [0,1), and r* as in ([4.11\) . 
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The entire discussion on the monotonic rational cubic spline FIFs can be encapsulated in the 
following theorem: 


Theorem 4.1. For a given set of monotonic data {{xi, i/i) : i = 1,2,N}, let S be the rational 
cubic spline FIF described in d3.3l) . Assume that the necessary conditions on the derivative 
parameters are satisfied. Then, the following conditions on the scaling factors and the shape 
parameters on each subinterval are sufficient for S to be monotonic on I = [xi,xn] '■ 


0 < a* < min 



d\ d]\[ 


* ^ihi 
K - 

vn - yi 


; [ 0 , 1 ), 


^ ^ hfdj + dj+i) - afxN - xi){di + dN) 

hiAi - afyN - yi) 

In particular, ifri = l -\— ^ ^ s reduces to the rational 

_ - ctifjN - 2/1 ) 

quadratic spline FIF given in d4.i2D . and in this case the above mentioned conditions on ai alone 
are sufficient for the rational FIF to be monotonic. 


4.2 Convex Data 

We assume a strictly convex set of data so that: 


Ai < A 2 < ■ ■ ■ < Atv-i- (4.13) 

To have a convex interpolant S and to avoid the possibility of S having straight line segments, it 
is necessary that the derivatives at knot points satisfy 


di < Ai < d2 < A 2 < ■ ■ ■ < di < Ai < ■ ■ ■ < dN- ( 4 . 14 ) 

For a concave data set inequality will be reversed. Let s be the classical counterpart of S studied 
elaborately in [fT9]| . Since s may fail to have second derivative at knot points, s is not twice 
differentiable on the entire interval /. Hence, in contrast to the claim made in [[T9l . convexity 
of s cannot be derived from the result that reads: s(x) is convex if and only if s^‘^\x) > 0 
for all X E I. However, the following result from elementary calculus justify the procedure 
adapted in [fT9l : Suppose that f is piecewise with increasing slopes, i.e., there is a subdivision 
a = xq < xi < ■ ■ ■ < Xk = b of I = [a,b] such that (i) f is continuous on I (ii) f is of class 
on each subinterval {xi-i,Xi), i = 1,2,... ,k (Hi) f has one-sided derivative at xi, X 2 , ■ ■ ■, Xk-i 
satisfying f(xf) < f(^x^) for i = 1,2,... ,k — 1, then f is convex on I. 

Now turning our attention to the rational cubic spline FIF S, it is worth mentioning that due to 
the fractality, S may not be even piecewise C^. Consequently, we cannot adapt the convexity 
procedure for the classical cubic spline s by applying the result stated above. Instead, we use the 
following results: 

(i) A differentiable function of one variable is convex on an interval if and only if its derivative is 
monotonically increasing on that interval. 
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(ii) Let f be a continuous function on [a, h]. If for each x G (a, h) one of the one sided derivatives 
y(i)(^+) gy exists and is nonnegative (possibly +ooj, then f is monotonic increasing 

on [a, h]. 

Owing to these results, to establish the eonvexity of S' e C^(/), it is enough to show that 
or ■) exists and is nonnegative for eaeh x G (xi, xn)- It is to this that we now turn. 

Informally, 


S^‘^\Li{x)) = ^S^^^(x) + Rf\x), 


(4.15) 


where e!^^\x) = + -6) + Q0(1 - Of + Dt{l - 0)^] 


hi[l + in - 3)6(1 - e)f 


A* = ri{di+i - Ai) + di- di+i - ^{rfd^ixM - a;i) - (vn - Vi)] + (xn - Xi)(di - d^)},^ 
B* = 3(di+i - Ai) - ^[d^ixM - xx) - (vn - Vi)], 

hi 

Q 

C* = 3{Ai - di) - -^[{Vn - yi) - di(xN - Xi)], 

D* = rfAi - di) + di- dj+i - ^{rfivN - Vi) - di(xN - xi)] + (xat - xi)(di - dTv)}. 

hi > 

We assume that 0 < «* < naj for i e J, where 0 < k < 1. Using the faet that Lj : [xi, x^v] —)■ 
[xj,Xj_|_i] satisfies Lj(xi) = Xj, Lj(xN) = x^+i, for j G J we get: 




a 


1 


+ ^; J = 2,3,..., iV - 1, 


9 4* 


iAT-l 




"Af-l 


(4.16) 


From (14.161) . it follows that the seeond derivative (right-hand) at the knot points Xj, j G J, and 
the seeond derivative (left-hand) at the extreme end point xjv is nonnegative if: D*,j G J, and 
Alj_i are nonnegative. For a typieal knot point Xj,j G J : 


S(2)(L,(x,)+) = ^S(2)(a;+) + /?f)(x+) 


(4.17) 


Assuming D*, j G J, to be nonnegative, (14.171) suggests that S*^^^ (Li(xj)+) > 0 is satisfied, 
provided i7.^^(x+) > 0. Again, i?®(x)^) > 0 is satisfied if: 


A* > 0, B* >0,C* > 0, and D* > 0; i G J. 

From the Three Chords Lemma for eonvex funetions [|5]|, it follows that a eonvex set of data 
should necessarily satisfy di < < ^at, where inequalities remain strict for strict convexity. 
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Now 5* > 0 di+i — Aj > ^[dN{xN — xi) — {yN — yi)]- Observing that if dN{xN — — 

(yN — yi) = 0, then B* > 0 is obviously satisfied, we get the condition on the scaling factor as 

^ Aj) 

* “ dN{xN - xi) - {yN - yi)' 

Similarly, C* >0 Ai- di> ^^[{yN - yi) - di{xN - a^i)] gives 


^ ^ ^i(Aj dj) 

* “ {yN - yi) - di{xN - xi)' 

Therefore, to obtain {Li{xj)~^) > 0 for alH G J and knot points Xj,j G J, it suffices to have 
A* > 0, D* > 0, and 


0 < a* < min < , 


di{di^i Aj) 


hi{Ai di^ 


dN{xN - Xi) - {yN - yi) ’ {yN - yi) - di{xN -Xi))' 
The following conditions on the shape parameter give A* > 0,D*>0. 

di+i - di + {ai/hi){xjsi - x{){di - djv) di+i - di + (ai/hi){x]^ - xi){di - djv) 


(4.18) 


Vi > max 


r ai+i — ai 
^ doJ-^ — Ao — I 


. di+i - Ai - (ai//ii)[djv(a:jv - xi) - (y^ - yi)\ ’ Ai - di - («i//ii)[(i/jv - yi) - di{xN - xr) 

The condition on stated above is equivalent to 

Mi 


Ti > 1 + 


rrii 


(4.19) 


where 

(y.' ol’ 

Mi = max{cii+i - Aj - -j^[dN{xN - xi) - {yN -yi)], di --j^[{yN-yi) - di{xN - xi)]}, 

rrii = min{cii+i - Ai - -j^[dN{xN - xi) - {yN -yi)], Ai-di - -j^[{yN - yi) - di{xN - xi)]}. 

Therefore the conditions (14.181) on the scaling factors and (14.191) on the shape parameters ensure 
(Li(a:j)+) > 0 for alH, j G J and S'(^)(a:“) > 0. Since the rational fractal function is 
generated recursively and [xi, x^v] is the attractor of the IPS {i?; Li{x) : z G J} (Li(xj)+) > 
0 for all i,j G J yield S'^^)(x+) > 0 for all x G (xi,Xiv). Hence, by the result quoted at 
the beginning of this section is monotonically increasing, as a consequence of which S is 
convex. Analogous procedure applies to a concave data set. 

Remark 4.5. If Ai = Ai_|_i, then for a convex fractal interpolant, we choose the scaling factor 
ai to be zero. Also, di = dj+i = A*. Thus, in this case the rational cubic FIF becomes S{x) = 
{xi+i-x)yi+{x-xi)yi+i ^ ^ ^ ^ ^ to a straight line segment on bj, Xj+i], as would 

^i+l L I J 

be expected. 
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Remark 4.6. When all a* = 0, the condition M.lSi is obviously true and the condition d4.i9l) 
reduces to 


M* 

r* > 1 H- with 

m* 


(4.20) 


M* = max{(ij+i - Aj, A* - di}, m* = min{(ij+i - Aj, A^ - di). 


This provides sufficient conditions for the convexity of the classical rational cubic spline [ 
p. 971 ]. 


Remark 4.7. In particular, if we choose 

Mj rrij 

ri = l-\ -h —, 

m,- Mi 


(4.21) 


with ai satisfying 0 < a* < min |a? ^ - hi{Ai di) - settle the 

convexity in question, then the rational cubic spline FIF S in AS.SD reduces to a lower-order form 
given by 

S(Li(x)) = aiS(x) + (1 - h)(pi - aiPi) + h(pi+i - aipjv) - 

(y.' 

where Gi := (dj+i - A*) - ■^[(iAr(xAr - Xi) - {ijn - Vi)] and 

CX' 

Hi := (Aj - di) - yUvn - yi) - di{xN - a;i)]. 

The classical counterpart of d4.22l) obtained by choosing all the scaling factors to be zero is 
described in / f7^ . In other words, d4.22D yields a fractal generalization of the classical rational 
spline with quadratic numerator and linear denominator studied in fiWl . 

Our particular choice of shape parameters given in d4.2iD verifying the convexity condition can 
be justified as follows. For the shape parameters as in M.2H . and the scaling factors satisfying 
|Q!j| < af, we have r* — 3 = 0{hf ), and consequently we obtain optimal 0{lT^) bound on 
interpolation error provided derivatives are estimated with 0{hf) accuracy. 


The main points in the above diseussion are extracted in the form of following theorem: 


Theorem 4.2. Given a convex (concave) data {(xi,yi) : i = 1,2,..., A}, assume that the 
derivative parameters satisfy the necessary convexity (concavity) condition. Then, the following 
conditions on the scaling factors and the shape parameters are sufficient for the corresponding 
C^-rational cubic spline FIF S to be convex (concave) on I = [xi, xtv]- 


0 < Oi < min l^a?, k 


Aj) 


1 




-|;k, G [0,1), 


n > max • 


dN{xN - Xi) - {VN - yi) ’ iVN - yi) - di{xN - Xi). 
di+i - di + (ai/hi){xN - xi){di - dN) dj+i - dj + {ai/hi)(xN - xi){di - ^at) 


di+i - Ai - {ai/hi)[dNixN - a^i) - (yN - yi)] ’ Aj - d* - (ai/hi)[yN - yi - di{xN - xi)] 
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4.3 Convex and Monotonic Data 

We now consider the possibility that the data satisfy both the monotonic increasing condition 
Vi < y 2 < ■ ■ ■ < Vn, and the strictly convex condition (14.131) . The derivative parameters must 
then satisfy the following inequalities: 


0 < di < Ai < ^2 < A 2 < ■ ■ ■ < Ai_i < di < Ai < ■ ■ ■ < (In- 


(4.23) 


We claim that the convex interpolation method described in the previous subsection is suitable 
for obtaining a convex and mono tonic fractal interpolant. To verify this claim, we proceed as 
follows. Assume that the sufficient conditions (14.181) on the scaling factors that achieve the 
convexity of the rational cubic FIF hold. Rearrangement of these inequalities gives 


~ dN{xN - Xi) - {un - yi) 


~ - yi) < di+l - ^dN{xN - Xi) 

lii hi 


and 


Q ^ di) 

^ {yj^ - yi) - di{xN - Xi) 


di - ^dii^N - Xi) < Ai 

rii 


hi 


(yN-yi)- 


Combining two inequalities obtained above, we have 

di - -^di{xN - xi) < Ai - -r^{yN - yi) < di+i - -^dN{xN - xi) (4.24) 

hi hi hi 

Since a* < 1, the condition Oj < /taf given in (14.181) implies the condition ai < Ktti in (14.91) . 
Also 


11 ' CX' 

ai < Kai,di > OVi G J aidi < k - —di di — -Adi{x]\f — Xi) > di — ndi (4.25) 

Xff — Xi hi 

Hence from (14.231) and (14.251) . we have di — ^di{xN — Xi) > 0. Consequently, (14.241) yield 
Ai — ^{yN — 2 / 1 ) > 0 and di+i — ^d^ixM — xi) > 0. Thus, we get the sufficient condition on 
the scale factors «*(? = 1,2 ,..., A — 1) that retain the data monotonicity. 

Assume that the sufficient condition (14.191) on the shape parameters r* for the convexity of the 
rational cubic FIF is true. We will prove that, this condition implies the condition (14.61) on r* for 
the monotonicity of the rational cubic fractal interpolant. Without loss of generality, assume that 

Mi = [di+i - -j^dN^XN - Xi)] - [Aj - -j^iyN - yi)], 
rrii = [Ai - ^{yN - yi)] - [di - ^di{xN - a:i)]. 

hi hi 

Denot e P* = di - f^di{xN - xi), Q* = A* - f^{yN - yi), R* = dj+i - f^dN^XN - xi). From 
(14.241) . we have P* < Q* < R*. Again, with these notations 


M, = R* - Q* = max{A: - QIQ* - P*} 


Q*< 


P* + R* 


(4.26) 
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The sufficient condition (14.61) for the monotonicity of a rational cubic FIF can be rearranged as 


> 


di + di+i-f:{di + dN){xN-xi) P* + R* 


_ hj, 

A.- 


hi 


Hvn -yi) 


Q*i 


(4.27) 


The sufficient condition (14.191) for the convexity of a rational cubic FIF in above notations be¬ 
comes 


ri>l + 


m-Q* 

Q* - P* 


R* - P* 
Ql-Pf 


(4.28) 


Note that (14.281) implies (14.271) . if which is equivalent to the condition described 

in (14.261) . But the condition (14.261) is obviously true due to our assumptions. The proof is similar 
if we assume that Mi = Q* — P* and rui = R* — Q*. Thus, we have proved the sufficient condi¬ 
tion for the convexity of a rational cubic FIF on shape parameters r j gives the sufficient condition 
on for the monotonicity r*. Therefore we conclude that for a given monotonic increasing con¬ 
vex data set, if derivative parameters are chosen according to (14.231) . then convex interpolation 
scheme developed in Section 4.2 will automatically produce a convex monotone rational cubic 
spline FIF. 


Theorem 4.3. Givenasetof strictly convex monotonic increasing data {{xi, Hi) : i = 1,2,. 
assume that the derivative parameters satisfy the necessary condition expressed in (14.231) . Then 
a convex interpolant obtained through the convexity preserving rational FIF scheme in Theorem 
(14.21) will automatically render a convex and monotone fractal interpolation curve. 


5 Some Remarks and Possible Extensions 

(i) Preserving Positivity. The proposed rational FIF can also generate positive fractal curves 

for a given set of positive data {(x*, ?/*) : i = 1, 2,..., A^}. Recall that the rational cubic 
spline FIF S is generated iteratively using the functional equation Si^Lf.)) = aiS{.) + 
Rf.). Hence with a* > 0 for all i E J, the conditions Rfx) > 0 for all x G / 
and for all z G J is enough to ensure S{x) > 0 for all x G /. Since Ri has pos¬ 
itive denominator, the positivity of Rfx) reduces to the positivity of cubic polynomial 
Pfx). Computationally efficient sufficient conditions for the positivity of Pfx) = P*{0) 
is given by At > 0, Bi > 0, Ci > 0, and A > 0. With simple calculations we ob¬ 
tain the following conditions on the IFS parameters: 0 < a, < min joj, and 

Tj > maxi—1, hidi+i-aidN(xff-xi) ^ particular, for all Oj = 0, we ob- 

tain conditions for the positivity of the rational cubic spline introduced in 11191 . It seems 
that lfT9l does not address the possibility of preserving positivity with the rational cubic 
spline developed therein. 

(ii) Admissibility of negative scalings for shape preserving: By taking monotonicity as an ex¬ 
ample of shape, we illustrate that the nonnegativity assumption on the scaling parameters 
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is not essential in our shape preserving schemes. For this purpose, we outline a slightly 
general problem, namely, identifying the parameters of the IFS so that the graph of 
lies in a prescribed rectangle R = I x [0,M]. Recall that 5(1) 

is generated using the 

IFS {(^Li{x), Fi^i{x,y))] where Fi^i{x,y) = By properties of the attrac¬ 

tor of the IFS, for 0 < S^^\x) < M it suffices to prove that Fj i(a;, y) G [0, M] for any 
(x, y) e R. Consider the two cases: (i) 0 < a* < a* (ii) —ai < ctj < 0. With 0 < a* < ai, 
0 < 1 / < M, the condition Fi^i{x,y) G [0,M] holds if 0 < Rf\x) < M{ai — a*). 
Now by the substitution of the rational expression rI^\x), the above inequality can be 
transformed to the positivity of suitable quartic polynomials. This provides conditions on 
IFS parameters for case (i). For case (ii), the condition Fi^i{x,y) G [0, M] will hold if 
—aiM < rI^\x) < ttiM. Again, this can be transformed to the positivity of suitable 
quartic polynomials. From this the conditions on the IFS parameters for case (ii) are de¬ 
duced. Combining the conditions in both cases, we obtain sufficient conditions on the IFS 
parameters for the graph of to lie in R. Taking M to be a large positive number, we 
deduce conditions for the monotonicity of S. Note that this allows negative values for the 
scaling parameters. 

(iii) Co-monotone/co-convex fractal interpolants: Often a data set will not be globally mono¬ 
tone, but instead switches back and forth between monotone increasing and monotone de¬ 
creasing. We need the interpolant to follow the shape of the data in the following sense: 
S is monotonically increasing/decreasing on a;*] if the data are monotonically increas¬ 
ing/decreasing on [xr,Xs]. Similarly we can define co-convex interpolation problem. If 
co-monotone/co-convex interpolation with the present fractal scheme is desired, a prelim¬ 
inary subdivision of the points into subsets of uniform shape is needed. Let us illustrate 
this with an example. Consider a data set {{xi, yi) : i = 1,2,... ,7} where yi < y 2 < f/s; 
1/3 > 1/4 > 1/5 > 1 / 6 ; He ^ yr- In order to achieve co-monotonicity, we must insure 
that slopes at transition points are zero, i.e., ds = 0 and cIq = 0. We divide the interval 
I = [xi,Xt\ into three subintervals such that the data possess same type of monotonicity 
property throughout that subinterval; R = [xi,X 3 \, R = [x 3 ,xe\, R = [xq,X 7 ]. We can 
apply the developed monotonicity preserving FIF algorithm to obtain a monotonically in¬ 
creasing rational cubic spline FIF Si on R. With proper renaming of the data points if 
necessary, the monotonically decreasing FIF algorithm can be applied to obtain a rational 
cubic spline FIF S 2 on R. Now consider the interval R = [xq, xf which contains only two 
data points. Here the iterations of IFS code cannot produce any new points. To remedy this 
problem, we introduce a new node say, (xl^ye) that is consistent with the shape present 
in [xq,Xj], i.e., Xq < xl < Xj and y% < yl < 1 / 7 . The “best” choice of additional node 
deserves further research. We apply the developed monotonically increasing FIF scheme 
with an arbitrary but shape consistent extra node to obtain a cubic spline FIF S' 3 , which is 
co-monotone with the data in R. Construct a rational cubic spline FIF S' in a piecewise 
manner by defining S|/. = Si. Note that the Hermite interpolation conditions on Si provide 
the Ckcontinuity for S. 

(iv) Optimal choice of parameters: Without a doubt, the scaling parameters and the shape pa- 
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rameters together provide a large flexibility in the ehoiee of an interpolant. Consequently, 
a natural question on an “optimal” ehoiee of the parameter arises. In this regard, some re¬ 
marks are in order. For higher irregularity in the derivative (quantified in terms of the fraetal 
dimension) we have to ehoose large values of the sealing faetors. In eontrast, for loealness 
of the seheme small values of ai are preferred. A preferable ehoiee of the sealing faetors 
and the shape parameters for a monotonieity/eonvexity preserving rational eubie spline FIF 
in terms of optimal interpolation error is given in Remarks (I4.4I) . (I4.7I) . Among the various 
shape preserving rational fraetal interpolants a visually improved solution may be obtained 
by minimization of a fairness funetional sueh as Holladay funetional. The feasible domain 
is given by suitable restrietion on the IFS parameters. This results in a eonstrained nonlin¬ 
ear optimization problem. In the classieal non-recursive shape preserving schemes, widely 

used Holladay functional is J (x)] dx. However, for the present fractal scheme, 

Xl 

may have discontinuity at each point of the interval / = [xi, xn], and consequently com¬ 
puting the integral occurring in this Holladay functional would be impossible at least in the 
Riemann sense. If we interpret the second derivative occurring in the Holladay functional 
loosely as (a;+) or assume that the FIF is in C^, then 


ieJ 


M = I ^ j + ^R^\Li\x)) dx, 


1 


h 




• 7 R 

l&J * 


+ Rq) 


where Rq = f R*{x) dx and R*{x) = \l~^{x)), if x G R. Applying the change of 

I 

variable x = L“Vx), 


M = YRm + R„^M= 

i-p,; 

i£j 


Thus, the constrained optimization problem is to Minimize M where variables are restricted 
according to finite set of inequalities resulting from the shape preserving constraints. It is 
felt that this constrained optimization problem can be solved by means of a differential 
evolution optimization algorithm/genetic algorithm. This procedure is justified if we make 
the interpolant to be by imposing suitable conditions on the derivative values and the IFS 
parameters resulting from the C^-continuity conditions. This may be done on lines similar 
to the cubic spline FIFs (see [I7]|, ifTTl l. This will lead to the fractal generalization of the 
standard C^-rational cubic spline introduced by Gregory [l24ll . 

From the point of view of approximation theory, the problem of finding optimal rational 
spline FIF S is an inverse problem which reads as: Given a set of values of a function, 
recover the IFS parameters generating this target function. Levkovich ff30ll has obtained 
contraction affine mappings generating a given function based on the connection between 
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the maxima skeleton of wavelet transform of the funetion and positions of the fixed points of 
the affine mappings in question. It is not without interest to note that for adapting a similar 
teehnique, the eonneetion between the strongest singularities of the FIF or its derivative and 
fixed points of the generating mappings, whieh is non- affine in the present ease, is to be 
developed rigourously. Lutton et al. ifSTl have applied genetie algorithms for solving this 
type of inverse problems. Resolution of the inverse problem is a major ehallenge of both 
theoretieal and practieal interest, whieh is not settled in its full generality. 

(v) A comparison of the present method with subdivision methods'. A possible alternative to 
the present reeursive shape preserving interpolant sehemes that introduee fraetality in the 
derivatives is so-ealled shape preserving subdivision sehemes (see, for instanee, [fT4ll20l 
[3^). Now we briefly eompare and eontrast the two methodologies. In both these interpola¬ 
tion methods, the desired interpolant is obtained eonstruetively. Convergenee of the fraetal 
interpolation scheme and the differentiability of the limit function follow from a straight 
forward application of the Banach contraction principle on a suitable function space. Es¬ 
tablishing the convergence of the scheme and the differentiability of the limit function is 
relatively harder in the subdivision schemes. Though the subdivision schemes add fraetality 
to the derivative function, we cannot directly control this fraetality in terms of the param¬ 
eters involved in the scheme. On the other hand, it is known |[T1 that as magnitudes of the 
scaling factors are increased from zero, the dimension of the derivative of the fractal spline 
increases. By controlling the scaling factors, the fraetality can be considered in a small 
portion of the domain, if in this part possible signal displays some complex disturbance. A 
quantitative measure of the irregularity (fraetality), namely, box counting/Hausdorff dimen¬ 
sion of the fractal curves in terms of the scaling parameters involved in the IFS is obtained 
in lUlIISl. Up to our knowledge, a quantitative measure of the irregularity of the deriva¬ 
tive in terms of the parameters involved is unavailable in subdivision schemes. Using the 
notions of hidden variable FIFs and coalescence hidden variable FIFs the present scheme 
can be extended to preserve shape of the data generated from a self-affine, non-self-affine, 
or partly self-affine and partly non-self-affine function. However, subdivision schemes do 
not specify about these properties of the constructed interpolants. The main appeal to the 
subdivision schemes resides in their localness. Due to the recursive and implicit nature of 
the FIF, the proposed scheme is, in general, non-local. However, the completely local clas¬ 
sical non-recursive interpolation scheme emerges as a special case of the proposed fractal 
interpolation method, and consequently, our method is local or global depending on the 
magnitude of the scaling factor in each subinterval. 


6 Numerical Examples 

Iterating the functional equation (cf. (13.31) ) with suitable choices of the scaling factors and the 
shape parameters as prescribed in Section S we generate different shape preserving rational cu¬ 
bic spline FIFs in this section. 

If the derivatives dj (1, 2,..., A^) are not supplied, estimates of derivatives are necessary. Meth- 
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ods that associate derivatives with data points involve estimates based on nearby slopes or data 
differences. Depending on the applications, various schemes based on linear combination (e.g., 
arithmetic mean method) or multiplicative combination (e.g., geometric mean method) of chord- 
slopes are developed in the literature (see, for instance, (IbUTTll ). With the notation Aj = , i e 

J, the three point difference approximation for the arithmetic mean method is given by di = 


hi—\-\-hi 


i = 2, 3,..., A — 1 with end conditions di = (1 -f ) Ai — 


hi A 


^ 3,1 — 


dN - yl + AtV-I - ^^An,N-2, An,N-2 


ys-yi 

^3 


yN-yN-2 
Xn—XN-2 ' 


6.1 Monotonicity Preserving Rational FIFs 

Consider a monotonic data set {{xi, Ui, di) = (0, 0,1.3333), (2,4, 2.6666), (3, 7, 2.6190), (9, 9, 
1.5833), (11,13, 2.4166)}. For monotonic FIFs, the computed bounds on the scaling factors are: 
0 < ai < 0.1818, 0 < 02 < 0.0985, 0 < 03 < 0.1538, 0 < 04 < 0.1818. 

Since S is obtained by iterating the functional equation (13.31) . perturbation of a particular scaling 
factor ttj and/or shape parameter may ripple through the entire configuration, i.e., interpolant 
is potentially non-local. However, we observed that the portions of the interpolating curve per¬ 
taining to other subintervals are not extremely sensitive towards changes in the parameters of a 
particular subinterval. To illustrate this, we take the monotonic rational cubic spline FIF in Fig. 
1 (a) as a reference curve, and analyze the effect of perturbing the parameters of a particular por¬ 
tion of this curve. Values of the parameters (rounded off to four decimal places) corresponding 
to various curves that are calculated according to the prescription in Theorem 14.11 are given in 
Tabled Changing the scaling parameter ai to 0.05 (see Table [B, we obtain Fig. 1(b). It is 
observed that the perturbation in ai effects the rational fractal interpolant considerably in the 
interval [xi, X 2 ], and there are no perceptible changes in other subintervals. Similarly, Fig. 1(c) 
and Fig. 1(d) are obtained by changing the scaling factor 02 and the shape parameter with 
respect to the reference curve. Effects of these changes are observed to be local. By taking zero 
scalings in each subinterval, we recapture a standard rational cubic spline due to Delbourgo and 
Gregory [|T9l (see Remarks [3Al 14.21) in Fig. 1(e). Optimal choices of the shape parameters sug¬ 
gested by Remark 14741 are used to generate the rational quadratic FIF in Fig. 1(f). 

Let us denote the monotonic rational cubic spline FIFs in Figs. I (a)-1(f) hy Si, i = 1, 2,..., 6 . 
Using the functional equation (14.11) . the derivative functions (i = 1, 2,..., 6 ) are generated 
in Figs. 2(a)-2(f). These curves possess varying irregularity. The derivative of the classi¬ 
cal rational cubic spline is smooth where as is nowhere differentiable. Note that has 
smoothness in the subinterval [x 2 ,xz] = [2,3] where the scaling factor is chosen to be 0. In 
this way, the fractality of can be restricted in a portion of the domain. It can be noted that 
due to the small values of the scaling factors in each interval is almost smooth. The frac¬ 
tal dimension of constitutes a numerical characterization of the geometry of the signal and 
may be used as an index for measuring the complexity of the underlying phenomenon (see, for 
instance, If36l l. 
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Table 1: Parameters corresponding to the Monotonic rational cubic FIFs 


Figure No 

Choice of parameters 

Fig. 1(a) 

di 

0.18 

0.09 

0.1 

0.18 

Ti 

2 

1.8 

31 

0.5 

Fig. 1(b) 

di 

0.05 

0.09 

0.1 

0.18 

Ti 

2 

1.8 

31 

0.5 

Fig. 1(c) 

di 

0.18 

0 

0.1 

0.18 

n 

2 

1.8 

31 

0.5 

Fig. 1(d) 

di 

0.18 

0.09 

0.1 

0.18 

Ti 

2 

1.8 

350 

0.5 

Fig. 1(e) 

Oti 

0 

0 

0 

0 

ri 

2 

1.8 

31 

0.5 

Fig. 1(f) 

Oii 

0.001 

0 

0.08 

0.001 

Ti 

2.9961 

2.7619 

23.8269 

2.9961 



(a): Monotonic rational cubic FIF 



(d): Monotonic rational cubic FIF 54 
(effect of rs) 



(b): Monotonic rational cubic FIF S2 
(effect of CKi) 



(e): Classical Monotonic rational cubic 
spline 55 


(c): Monotonic rational cubic FIF Sz 
(effect of 0.2) 



(f): Monotonic rational quadratic FIF Sq 


Figure 1: Monotonic rational cubic/quadratic spline FIFs: S'i, i = 1, 2,..., 6. 
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(a): Fractal function 5"^^^ 
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(c): Fractal function (smooth in 
[2,3]) 



(d): Fractal function (e): Function (f): Function 

Figure 2: Derivatives of monotonic rational FIFs in Figs. l(a)-l(f) 


6.2 Convexity Preserving Rational FIFs 

Consider the convex data set {(2.2, 2), (4, 0.625), (5, 0.4), (10,1), (10.22,1.8)}. We computa¬ 
tionally generate convex rational cubic spline FIFs by using (13.31) and the parameter values given 
by Theorem 14.21 The derivative parameters required for the implementation of the IFS scheme 
are estimated using the arithmetic mean method. The convex rational cubic spline FIF gener¬ 
ated in Fig. 3(a) is taken as a reference curve. Changing the scaling factor as to 0.005 and 
keeping the values of the other parameters as in Fig. 3(a) (see Table |2l), we obtain the convex 
rational cubic spline FIF in Fig. 3(b). It can be observed that the change in as influence the 
curve only in Further, due to a small value of the scaling factor and a large value of 

the shape parameter the FIF converges to a line segment in [xs, X 4 ], demonstrating the tension 
effect. Similar experiments may be conducted by changing the scalings in other subintervals 
and the shape parameters. By taking all the scaling factors to be zero and the shape parameters 
according to (14.211) . a classical rational quadratic spline that retain the data convexity is obtained 
in Fig. 3(c). Thus, Fig. 3(c) provides a numerical example for the convex rational quadratic 
spline by Delbourgo [fT^ . As in the monotonicity case, it can be observed by plotting the graph 
of the derivatives that the scaling factors provide fractality in the derivatives or (more 
precisely right hand second derivative). 


26 













































































Table 2: Parameters eorresponding to the eonvex rational cubic FIFs 


Figure No 

Choice of parameters 

Fig. 3(a) 

ai 

ri 

0.02 0.001 0.16 0.007 

37 8 269 8 

Fig. 3(b) 

CXi 

Ti 

0.02 0.001 0.005 0.007 

37 8 269 8 

Fig. 3(c) 

Oti 

Ti 

0 0 0 0 

3 2.6459 12.8069 3 





(a): Convex rational cubic FIF 


(b): Convex rational cubic FIF (effect of (c): Classical convex rational quadratic 


as) 


spline 


Figure 3: Convex rational cubic/quadratic spline FIFs. 


6.3 FIFs with Mixed Shape Properties 

The theoretical discussion we had in section!^ was confined to data with same shape character¬ 
istics in the entire interpolation interval. However, we can also apply our schemes with proper 
modification and mixing to obtain fractal interpolants for data with mixed shape properties. We 
illustrate this with two examples. 

The first example taken from lf32]| is a data set generated from a function $1 defined on [0,4] 
which is positive on [0,1], strictly increasing on [1, 2], constant on [2, 3], and concave on [3,4]. 
Suppose that we want to use the proposed fractal interpolation scheme to construct a approxi- 
mantto this function interpolating the data set{(0,1), (0.5, 0.2), (1,1), (1.5,1.3), (2,1.8), (2.5,1.8), 
(3,1.8), (3.5,1.5), (4, 0)} with the same shape characteristics. To achieve this, derivative values 
that are consistent with the required shapes are estimated: di = —3.2, d 2 = = de = df = 0, 

ds = 1.1, di = 0.8, ds = —1.8, and dg = —4.2. We divide the interval into four subintervals of 
unit length such that in each of the subinterval data possess the same shape property. To obtain 
a positive rational cubic spline FIF Si in [0,1], we iterate the functional equation (13.31) with the 
scaling factors and the shape parameters satisfying the required conditions (see Section |5] (i)). 
Our specific choices of the scaling parameters and the shape parameters are: ai = a 2 = 0.15; 
ri = 1.5, r 2 = 0.5. On the interval [1, 2], we apply our monotonicity preserving algorithm with 
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the parameter values ol\ = a2 = 0.3, ri = 8, and r 2 = 1 to obtain the monotonic rational cubic 
spline FIF S' 2 . Note that here the parameters are indexed by considering the interpolation to take 
place in the subinterval [1,2], not the entire interval. Following Remark (14.11) . we generate a 
linear interpolant S'3 on [2, 3]. Finally, the concavity preserving rational cubic spline FIF 5*4 is 
obtained on [3,4] by iterating functional equation (13.31) with parameters values satisfying con¬ 
cavity condition, our specific choices being ai = 0.15, 0^2 = 0.2, ri = 9, and r 4 = 4. Since 
the data satisfy both the monotonic decreasing condition and the strictly concave condition on 
[3,4], and derivative parameters are selected accordingly, the concave interpolation scheme au¬ 
tomatically render a concave and monotonically decreasing interpolant (see Section |431) . The 
fractal function S on [1,4] (see Fig. 4(a)) is obtained by pasting Si, i = 1, 2, 3,4. Since Si and 
S^^'^ (i = 1, 2, 3,4) are continuous, the continuity of S and follows from the pasting lemma. 
Consequently, the fractal function S G C^[0,4] given in Fig. 4(a) provides an approximation to 
$1 satisfying the required shape properties. 

Consider a function ^2 defined on [0,10], which is monotonic decreasing on [0,4], monotonic 
increasing on [4, 6], and monotonic decreasing on [6,10]. Our second example is concerned with 
the construction of a -approximation which is co-monotone with this function, where all we 
know about the function is the function values at specified points, say, {(0,10), (1.5, 5), (4, 3.5), 
(6, 7.1), (8, 3), (10, 0)}. As discussed in section|5l we divide the interval in to three subintervals 
Ji = [0,4], I 2 = [4,6], and Ja = [6,10] where the data are monotonic decreasing, monotonic 
increasing, and monotonic decreasing respectively. Since the interval I 2 = [4, 6] contains only 
two knot points, the FIF scheme demands insertion of a node in this interval. Let the new node be 
(5, 6). Derivative values that are consistent with the required shapes are chosen as di = —4.35, 
d 2 = —2.31, ds = 0, d^ = 1.8 (at the inserted knot), d^ = 0, d^ = —1.77, and dj = —1.2. A 
rational cubic spline FIF Si is constructed on Ji by taking ai = a 2 = 0.2, ri = 2, r 2 = 12, 
and iterating the functional equation (13.31) . On I 2 , the functional equation (13.31) with the param¬ 
eter values Q!i = Q !2 = 0.3, ri = 2, and r 2 = 91 generates S' 2 . Finally iterations of (13.31) with 
oi = «2 = 0.4, ri = 2, and r 2 = 26 yield S3. The fractal function S G defined in a piecewise 
manner by SL = Si, z = 1, 2, 3 is co-monotone with $ 2 , and it is given in Fig. 4(b). 

1 Conclusions 

A new kind of rational cubic fractal splines involving shape parameters is proposed in the present 
work to provide a tool for univariate shape preserving interpolation. Number of parameters in 
the rational IFS is kept to be minimum (one family of parameters for controlling fractality in 
the derivative function and one for providing shape preserving characteristics) for computational 
efficiency. Due to the presence of the scaling factors and the shape parameters involved in the 
definition, the proposed C^-rational cubic spine FIF generalizes the classical rational splines stud¬ 
ied in the references [I16U19U25II . Despite the implicit and recursive nature of the FIFs, it is shown 
that the existence of range restricted fractal interpolants depends only on the solvability of a fi¬ 
nite set of inequalities resulting from the constraints. These inequalities are shown to be solvable 
if the shape parameters are above and the scaling factors are below certain explicitly calculable 
bounds. Uniform convergence of the rational cubic spline FIF to the original data generating 
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(a): Rational cubic spline FIF with mixed shape property (b): Co-monotone rational cubic FIF 

Figure 4: Rational cubic spline FIFs with mixed shape properties. 


function / G C^{I) is established. The convergence analysis shows that 0{K'){r = 1,2, 3,4) 
error bounds can be achieved by suitable choices of the derivatives, the scaling factors, and the 
shape parameters. Thus, the present interpolation method has convergence properties similar to 
that of its classical counterpart, which should be considered along with the flexibility and diver¬ 
sity offered by the new method. The scaling factors and the shape parameters can be selected 
suitably to find an interpolant satisfying chosen properties such as smoothness, approximation 
order, locality, fractality in the derivative, and shape preservation of the data. The fairness (visual 
pleasantness) of the interpolant can be achieved through a constrained non-linear optimization. 
For the shape preserving interpolants with varying irregularity in the derivatives, the result is 
encouraging for the fractal spline class treated in this paper. Consequently, it is felt that the pro¬ 
posed scheme can provide an efficient mathematical tool for the simulation of curves occurring 
in the study of physical systems, for instance, in the study of nonlinear control problems such as 
pendulum-cart system and in some fluid dynamics problems such as motion of a falling sphere 
in a non-Newtonian fluid. 
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